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Most of the mathematical approaches for quantum 
Langevin equation are based on the non-commutativity of 
the random force operators. Non-commutative random force 
operators are introduced in order to guarantee that the equal- 
time commutation relation for the stochastic annihilation and 
creation operators preserves in time. If it is true, it means that 
the origin of dissipation is of quantum mechanical. However, 
physically, it is hard to believe it. By making use of the uni- 
fied canonical operator formalism for the system of the quan- 
tum stochastic differential equations within Non-Equilibrium 
Thermo Field Dynamics, it is shown that it is not true in 
general. 



I. INTRODUCTION 

The studies of the Langevin equation for quantum sys- 
tems were started in connection with the development 
of laser and are still continuing in order to de- 

velop a satisfactory formulation |^-^| (see comments in 
|jiof ) . Most of the mathematical approaches for quantum 
Langevin equation are based on the non-commutativity 
of the random force operators. For dissipative systems, 
for example, we have equations for the operators (a(t)} 
and {a?{t)) averaged with respect to random force oper- 
ators of the forms 



-<a(t)) = -iw<<»(*)> -*<<*(*)>, 



dt 



(a^t)}=itu{a\t)}- K (a^{t)}, 



with the initial condition 



(a(0))=a, (^(0))=^, 



(1) 
(2) 

(3) 



where a and a 1 satisfy the canonical commutation rela- 
tion 



[a, aJ] 



1. 



(4) 



The equal-time commutation relation for these operators 
decays in time: 



Random force operators df(t) and dp(t) are introduced 
in order to rescue this situation. If the random force 
operators in the Langevin equations 



da(t) = -iua(t)dt - na(t)dt + V2k df(t), (6) 
da\t) = iuja\t)dt - Ka){t)dt + \/2k dp(t), (7) 



satisfy 



W{t), df(t)] = dt 



(8) 



the equal-time commutation relation for the stochastic 
operators a(t) and a'(t) preserves in time: 



meaning that 



d([o(t), ot(t)]) =0, 



[a(t), o+(t)] = l, 



(9) 



(10) 
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with (|). 

The above argument is of zero temperature related 
only to the zero-point fluctuation. However, it has been 
extended to include the situations for finite temperature. 
Then, we have a crucial question. Should we interpret 
that the origin of thermal dissipation is quantum me- 
chanical? In this paper, we will investigate this question 
with the help of the system of the stochastic differential 
equations within Non-Equilibrium Thermo Field Dynam- 
ics (NETFD) 

NETFD is a canonical operator formalism of quantum 
systems in far-from-equilibrium state which provides us 
with a unified formulation for dissipative systems by the 
method similar to the usual quantum field theory that 
accommodates the concept of the dual structure in the 
interpretation of nature, i.e. in terms of the operator al- 
gebra and the representation space. The representation 
space of NETFD (named thermal space) is composed of 
the direct product of two Hilbert spaces, the one for non- 
tilde fields and the other for tilde fields. It can be said 
that NETFD is a framework which gives a foundation of 
Green's function formalisms, such as Schwinger's closed- 
time path method, Keldysh-method, and so on [^9| [fTf , 
in terms of dissipative quantum field operators within the 
representation space constructed on an unstable vacuum. 

In the extension to take account of the quantum 
stochastic processes p3|-p8fl , NETFD again allowed us to 
construct a unified canonical theory of quantum stochas- 
tic operators. The stochastic Liouville equations both of 
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the Ito and of the Stratonovich types were introduced in 
the Schrodinger representation. Whereas, the Langevin 
equations both of the Ito and of the Stratonovich types 
were constructed as the Heisenberg equation of motion 
with the help of the time-evolution generator of corre- 
sponding stochastic Liouville equations (Fig. 0) . The Ito 
formula was generalized for quantum systems. 

NETFD has been applied to various systems, e.g. 
the dynamical rearrangement of thermal vacuum in su- 
perconductor p9| , spin relaxation , various tran- 
sient phenomena in quantum optics ]31-|35|], non-linear 
damped harmonic oscillator (3(|, the tracks in the cloud 
chamber (a non-demolition continuous measurement) 
p7[ , microscopic derivation of the quantum Kramers 
equation |38|. 



section |jv|, two systems will be applied to the model of 
damped harmonic oscillator interacting with irrelevant 
random force system by the linear dissipative coupling. 
It will be shown that both systems are consistently ap- 
plicable. The existence of the non-commutative random 
force operators is essential for the system with unitary 
time-evolution generator. In section [v| two systems will 
be applied to the model of damped harmonic oscillator in- 
teracting with irrelevant system by the position-position 
interaction. It will be shown that the system with uni- 
tary time-evolution generator cannot produce the frame- 
work which is consistent with the master equation, since 
there appear only commutative random force operators 
in martingale. Section VI will be devoted to summary 
and discussion. 



II. FRAMEWORK OF NETFD 



Dissipative Q. Field Th. 



Ave. Eq. of Motion 



VE 



-A{t) = i[H(t), A(t)] 



VE 

0(t)) = -iH\0(t)) 



dA(t) = i[H f (t)dt, A{t)] 
-dM[t) [dM(t), A(t)] 



Langevin Eq. 
(Ito) 



Heisenberg Eq. 
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d\O f (t)) = -iH f ,tdt |0/(t)) 
Hf,tdt = Hdt + dM t 



Stoch. Liouville Eq. 
(Ito) 



Heisenberg Rep. 



Schrodinger Rep. 



FIG. 1. System of the Stochastic Differential Equations 
within Non-Equilibrium Thermo Field Dynamics. RA stands 
for the random average. VE stands for the vacuum expecta- 
tion. 

In the next section, the framework of NETFD will 
be briefly explained. In section [II, two systems of the 



stochastic differential equations will be introduced, one 
with non-unitary time-evolution generator and the other 
with unitary time-evolution generator. Both systems are 
constructed to be consistent with the same quantum mas- 
ter equation. The key is the existence of the fluctuation- 
dissipation theorem between the multiple of martingale 
operators and the imaginary part of hat-Hamiltonian. In 



The dynamics of physical systems is described, within 
NETFD, by the Schrodinger equation for the thermal 
ket-vacuum \0(t)): 



^\0(t)) = -iH\0(t)). 



(11) 



The time-evolution generator H is an tildian operator 
satisfying 



(iff)~ = iH. 

The tilde conjugation ~ is defined by 

(AiA,)"=U, 
(ci^i + c 2 A 2 y = c\A x + c* 2 A 2l 
(A)~ = A, 
(A*)- = it, 



(12) 

(13) 
(14) 
(15) 
(16) 



where c\ and c 2 are c-numbers. The tilde and non-tilde 
operators at an equal time are mutually commutative: 



[A, B] = 0. 



(17) 



The thermal bra- vacuum (1| is the eigen- vector of the 
hat-Hamiltonian H with zero eigen- value: 



a\H = o. 



(18) 



This guarantees the conservation of the inner product 
between the bra and ket vacuums in time: 

(l|0(i)) = 1. (19) 

Let us assume that the thermal vacuums satisfy 

(1|~ = <1|, |0(t )>~ = |0(*o)), (20) 



at a certain time t = to. Then, ( |12[ ) guarantees that they 
arc satisfied for all the time: 
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(II 



|Q(t))~ = |0(t)>. 



(21) 



The tilde operator and the non-tilde operator are re- 
lated by the thermal state condition for the bra vacuum: 



(i\A 



(22) 



which reduces the numbers of the degrees of freedom to 
the original ones. The numbers of the degrees of freedom 
were doubled by the introduction of tilde operators. 

The observable operator A should be an Hermitian op- 
erator consisting only of non-tilde operators. 



III. TWO SYSTEMS OF STOCHASTIC 
DIFFERENTIAL EQUATIONS 

A. Quantum Master Equation 

Let us consider the system of quantum stochastic dif- 
ferential equations which is constructed to be consistent 
with the quantum master equation (the quantum Fokker- 
Planck equation) 



d_ 

di 



\0(t)) = -iff|0(t)>, 



with the hat-Hamiltonian 

H = H s + in, 

where 

Hs = Hg - Hs , 



(23) 



(24) 



(25) 



with Hs being the Hamiltonian of a relevant system. It 
is easily seen that H$ satisfies 



(1\H S = 0. 



(26) 



It is assumed that the imaginary part 77 of the hat- 
Hamiltonian can be divided into two parts, i.e., the re- 
laxational part 77r and the diffusive part Tip: 



n = n R + n D , 

and each of them satisfies 

(l|77 fl = 0, (114b = 0. 
Introducing the time-evolution operator V(t) by 
d 



dt 



v(t) = -mv{t), 



(27) 
(28) 

(29) 



with the initial condition V(0) = 1, we can define the 
Heisenberg operator 

A(t) =V- 1 {t)AV{t), (30) 

which satisfies the Heisenberg equation 



-A(t)=i{H(t), A(t)], 



(31) 



for dissipative systems. 

The equation of motion for the averaged quantity 
(l|A(t)|0) is derived by means of the Heisenberg equa- 
tion (J3l|) by taking its vacuum expectation: 



f t (l\A(t)\0)=i(l\[H(t), A(t)]\0). 



(32) 



The same equation can be also derived with the help of 
the master equation (23) as 



-i(l\AH\0(t)). 



(33) 



We would like to emphasize here that the existence of 
the Heisenberg equation of motion ( |3l| ) for coarse grained 
operators is one of the notable features of NETFD. This 
enabled us to construct a canonical formalism of the dis- 
sipative quantum field theory, where the coarse grained 
operators a(t) etc. in the Heisenberg representation pre- 
serve the equal-time canonical commutation relation 



[a(t), t(f)] = 1, 



l(t), at(t)] = i 



(34) 



Note that we have an equation of motion for a vector 
(l\A(t): 



-(l\A(t)=i(l\[H(t), A(t)] 



i(l\[H s (t), A(t)} 
-n{(l\[A(t), tf(t)]a(t) 
+<l|ot(i)[o(t), A(t)}} 
+2Kfi(l\[a(t), [A(t), a\t)]] 



(35) 



in terms of only non-tilde operators with the help of the 
condition (fn]). Applying the ket- vacuum |0) to (|35|), we 
obtain the equation of motion for the averaged quantity 



B. Non-Unitary Time- Evolution 

The system of stochastic differential equations with 
non-unitary time-evolution is constructed by the follow- 
ing general procedures. 

The stochastic Liouville equation 



d\0 f (t))) = -iH f<t dt\0 f (t))), 



(36) 



of the Ito type is specified with the stochastic hat- 
Hamiltonian 



Hf, t dt = Hdt + dMt 
= H s dt + illdt 



dM f 



(37) 



where 77 = 77^ + 77£> is the same that appeared in the 
master equation (23). The martingale operator dM t an- 
nihilates the bra- vacuum (1| of the relevant system: 
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(l\dM t = 0, 



(38) 



Hfjdt = H s dt + iII R dt + dM t , 



(48) 



which means that the stochastic Liouville equation (|3q ) 
preserves its probability just within the relevant sys- 
tem. This feature is the same as the one within the sys- 
tem of stochastic differential equations for classical sys- 
tems. The martingale operator satisfies the fluctuation- 
dissipation theorem of the second kind: 



dM t dM t = -2ir D dt, 

which should be interpreted as a weak relation.]] 
Applying (| to (|36|), we have an equation for 

|0(*)) = (|0/(t)», 



(39) 



(40) 



which is nothing but the quantum master equation (23) 
Introducing the stochastic time-evolution operator 



dv f (t) = -m u dtVf(t), (4i) 

we can define the stochastic Heisenberg operator 

A(t)=Vf\t)AV f (t), (42) 

which satisfies the stochastic Heisenberg equation (the 
Langevin equation) 



dA{t) = i[H f (t)dt, A(t)} 

-dM(t) [dM(t), A(t)] 



(43) 



of the Ito type. Here, we introduced the martingale op- 
erator in the Heisenberg representation by 



Note that 



with 



dM{t) = d (Vj- x {t)M t Vf{t) 



dM(t) = d'M(t), 



d'M(t) = Vf^^dMtVjit). 



(44) 



(45) 



(46) 



Making use of the relation of the Ito-Stratonovich 
stochastic calculus (see Appendix ]b|) , we can derive from 
( |36| ) the stochastic Liouville equation 



d\0 f (t))) = -iHf, t dto\0 f (t))), 



(47) 



of the Stratonovich type, where the symbol o indicates 
the Stratonovich stochastic multiplication. The stochas- 
tic hat-Hamiltonian 



It is similar to the classical cases where the fluctuation- 
dissipation theorem of the second kind is specified within the 
stochastic limit. 



contains only the relaxational part 77^. 

With this hat-Hamiltonian, we can write down the 
stochastic Heisenberg equation 



dA(t) = i[H f (t)dt ° A(t)], 



(49) 



of the Stratonovich type. Note that it does not have 
the term producing diffusive time-evolution, which is the 
same characteristics that appeared in the system of clas- 
sical stochastic differential equations. 



C. Unitary Time- Evolution 

The system of stochastic differential equations with 
unitary time-evolution is constructed by the following 
general procedures. 

The stochastic Liouville equation 



d|0,(t)» = o|<v(t)», 



(50) 



of the Stratonovich type is specified with the stochastic 
hat-Hamiltonian 

H% t dt = H s dt + dMf, (51) 

with the Hcrmitian martingale operator 



(dM?y 



'dM t u \ =dM?. 



(52) 



Note that ( p0| ) does not satisfy generally the conservation 
of probability just within the relevant system, i.e., 



{\\dM? £ 0, 



(53) 



but it does within whole the system, the relevant and 
irrelevant systems, i.e., 



(l\dM t u = 0. 



(54) 



Here, = (1|(| with (| being the bra-vacuum of the 
quantum Brownian motion (see Appendix ^) . The mar- 
tingale operator satisfies the fluctuation-dissipation the- 
orem 



dM^dM^ = -2hdt, 



(55) 



of the second kind. 

Introducing the unitary stochastic time-evolution op- 
erator Uf(t) by 



dU f (t) = ~iHf t dtoU f (t), 



(56) 



with the initial condition Uf(0) — 1, we can define the 
stochastic Heisenberg operator 



A{t) = Uj l {t)AU f {t), 



(57) 
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which satisfies the stochastic Heisenberg equation (the 
Langevin equation) 



dA(t) = i[Hf(t)dt ° A(t)}, 



(58) 



of the Stratonovich type. Note that the time-evolution 
generator Uf(t) is a unitary operator: 



u}(t) = uj\t). 



(59) 



By making use of the relation between the Ito and 
Stratonovich stochastic calculus, we can derive from j50| ) 
the stochastic Liouville equation 

d\O f (t))) = -iH%db\O f (t))), (60) 

of the Ito type with the stochastic hat-Hamiltonian 

H% = Hdt + dM? '. (61) 



Applying (| to (60), we see easily that it reduces to the 
quantum master equation (|25|). 

Within the Ito calculus, the time-evolution operator 
Uf(t) satisfies 



dU f (t) = -iH^ it dtUf(t), 



(62) 



with the initial condition Uf(0) = 1. The stochastic 
Heisenberg operator A(t) defined by fl57| ) satisfies the 
stochastic Heisenberg equation 



dA{t)=i[H u f {t)dt, A(t)] 

-dM u {t) [dM u {t), Alt)], 



(63) 



of the Ito type. Here, we introduced the martingale op- 
erator in the Heisenberg representation by 



dM' 



Note that |37| 



with 



T {t)=d(uj\t)M?U f {t) 



dM u (t)=d'M u (t) 



d'M u (t) = Uj l {t)dM^U } {t). 



(64) 

(65) 
(66) 



IV. APPLICATION TO QUANTUM DAMPED 
HARMONIC OSCILLATOR 

A. Quantum Master Equation 



The hat-Hamiltonian of the semi-free field is bi-linear 
in (a,a,a' ,5J), and is invariant under the phase trans- 
formation a — > ae : 

H = g\a) a + <? 2 a^a + 930,d + g±a*(v + go, (67) 



where g's are time-dependent c- number complex func- 
tions. 

The operators a, oJ, etc. satisfy the canonical commu- 
tation relation]^] 



[a k , Oj c ,]=5k,k' ) [5k, a k /]=<5k,k'- 



(68) 



The tilde and non-tilde operators are mutually commu- 
tative. Throughout this paper, we do not label explicitly 
the operators a, a\ etc. with a subscript k for specifying 
a momentum and/or other degrees of freedom. However, 
remember that we are dealing with a dissipative quantum 
field. 

The tildian nature {iH)~ — iH makes ( p"7l ) tildian: 



H = u){a) a — a) a) + ill , 



(69) 



with 



II = c\{a)a + a) a) + c 2 aa + c^a 1 + 04, (70) 



tat 



where uj = 3?e g\ = — Sfte gi, c\ — 3m g\ — 3m g2, c 2 = 
3m (73, C3 = 3m 34 and C4 = 3m go. 
With the help of @ for A = a: 

(l|5=(l|at, (71) 

the property (1\H = gives us the relations 

2ci + c 2 + c 3 = 0, c 3 + c 4 = 0. (72) 

Then, (ffol) reduces to 



II = ci(a} a + a 1 a) 

—c 2 ad — (2ci + c 2 ) a) a) + (2ci + c 2 ) . (73) 

Let us write down here the Heisenberg equations for a 
and a): 



dt 



a(t) = -iua(t) + cta(t) - (2ci + c 2 ) a tt (t), (74) 



eft 



a tt (t) = icja tt (i) - cia tt (t) - c 2 a(i). 



(75) 



Since the semi-free hat-Hamiltonian H is not necessarily 
Hermite, we introduced the symbol jf in order to distin- 
guish it from the Hermite conjugation f . However in the 
following, we will use f instead of ff, for simplicity, unless 
it is confusing. By making use of the Heisenberg equa- 
tions ([74] ) and ( [75^ , we obtain the equation of motion for 
a vector (l|at(f)a(f) in the form 



d 
dt 



(l\ a i{t)a(t) = -2«(l|a 1 (i)a(t) + iS<{l\, (76) 



2 Throughout this paper, we confine ourselves to the case 
of boson fields, for simplicity. The extension to the case of 
fermion fields are rather straightforward. 
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where we introduced k and S < respectively by 



K 



ci + c 2 , 
i(2ci + c 2 ) 



(77) 
(78) 



In deriving (|76|), we used the thermal state condition ( |22] ) 
in order to eliminate tilde operators. 

Applying the thermal ket vacuum |0) to ([76|), we obtain 
the equation of motion for the one-particle distribution 
function 



as 



n(t) = (l|a tt (i)a(i)|0) = (l|at O |0(*)), 



—n(t) = -2nn{t) + iS < . 
at 



(79) 



(80) 



The equation ( p0[ ) is the Boltzmann equation of the sys- 
tem. The function S < is given when the interaction hat- 
Hamiltonian is specified. 

The initial ket-vacuum |0) = |0(i = 0)) is specified by 



a|0) = /ot|0>, 



(81) 



with a real quantity /. Here, we are neglecting the initial 
correlation p2| . The initial condition of the one-particle 
distribution function n = n(t = 0) can be derived by 
treating (l|aa|0) as follows. In the first place, 



(l|o5|0) = (l\afa?\0) 

= /((l| a t a | ) + (l|0)) 
= /(n + l), 



(82) 



where we used the tilde conjugate of ( pl|) for the first 
equality, and the canonical commutation relation ( |68| ) 
for the second. On the other hand, 



(l\aa\0) = (l\aa\0) 
= <l|o t a|0> 
= n. 



(83) 



Here, for the first equality, we used (|l7|), i.e., the commu- 
tativity between the tilde and non-tilde operators, and, 
for the second equality, (fn)). Equating ( [84 ) and (|83|), we 
see that 



1-/' 



/ = 



1 



(84) 



If it is assumed that there is only one stationary state, 
we can refer the stationary state as a thermal equilib- 
rium state. We will assign the thermal equilibrium state 
to be specified by the Planck distribution function with 
temperature T: 



n(t — > oo) = 
Then, we have from (j80|) 



1 



(85) 



iS < = 2nn. 
In this case, the Boltzmann equation 



— n(t) = -2/e(n(t) - n) 



(86) 



reduces to 



(87) 



substituting 



Solving ([77 ) and (78) with respect to c\ and C2, and 



(8tj) for S < into (|70|), we finally arrive at 
the most general form of the semi-free hat-Hamiltonian 
H corresponding to the stationary process p3[ : 



H = H, 



where 



H s = He 



H s , 



in, 



Hs = ujoja, 



H = -K [(1 + 2n) (a f a + a f a) 

-2 (1 + n)aa~ 2^^] - 2t 



(88) 
(89) 

(90) 



This hat-Hamiltonian is the same expression as the one 
derived by means of the principle of correspondence 
|IT| , ^2| when NETFD was constructed first by referring to 
the projection operator formalism of the damping th eory 
@||||. The hat-Hamiltonian (|f) with (||) and @) 
describes time-evolution of the system of a damped har- 
monic oscillator. 

A detailed investigation of the system is given in Ap- 
pendix^. Here we just write down an attractive expres- 
sion which leads us to a new concept, named spontaneous 
creation of dissipation. The expression is given by 



|0(*)> = cxp 



- / d 3 fc(l| 7M 7M|0) 7 M 



where 



(1|7M7M|0) = -n k (t)+n k (0) 



|0) (91) 



(92) 



is the order parameter for dissipative time-evolution of 
the unstable vacuum. The annihilation and creation op- 
erators 



u=l 

-u=l i 

%,t =7fc = 



H=2 ~l 

%t =7* S 

-n=2 

%,t = -7M) 



(93) 
(94) 



are defined by 

lit = B k (tra%, % = alB^ity^ (95) 
with the time- dependent Bogoliubov transformation: 



Bk[tr = fl + Mt) -nj(t) 



(96) 



They satisfy the canonical commutation relation 

Wm> r k ' <t } = K^, (97) 
and annihilate the bra- and ket-vacuums at time t: 



7mI0(*)> =0, 



(l|7f = 0- 
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B. Non-Unitary Time- Evolution 



d([a(t), t(t)]) =0. 



(Ill) 



Confining ourselves to the case where the interaction 
hat-Hamiltonian between the relevant system and the ir- 
relevant system of Brownian motion is bi- linear in (a, a' , 
and their tilde conjugates) and (dB t , dB\, and their tilde 
conjugates), and is invariant under the phase transfor- 
mation a — ► ae 10 , and dBt ^ dBt e 10 , the martingale 



operator satisfying (38) is given by 

: dM t : = i (i*dW t + J*dWt 



where we introduced the random force operator 
dW t 



^2k (jiABt + vdB[ 



and the annihilation and creation operator 



[ia + va' , 
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,t 



(99) 



(100) 



(101) 



of the relevant system with ji + v = 1, which satisfy the 
commutation relation 



and annihilate the relevant bra-vacuum: 



(1| 7 * = Q, (1|7^ 



0. 



(102) 



(103) 



Note that the normal ordering : • • • : in (]99) is defined 
with respect to the annihilation and the creation oper- 
ators. Making use of the annihilation and the creation 
operators, we can rewrite 77r and II o consisting of 77 
introduced in (|9^) as 

n R = -k (7 V + 7*>) , (104) 
II D = 2k (n + j/) 7*7*. (105) 
The Langevin equations for a(t) and a^(t) are given by 



da(t) — ~ioja{i)dt 

-k[(jjl - v)a(t) + 2va)(t)]dt + dW t , 
da f (i) = iuja){t)dt 

-K[2fj,a(t) - {p.- v)a)[t)]dt + dW u 

where we used the facts 

dW{t) = dW u dW(t) = dW t . 



(106) 
(107) 

(108) 



If we put v = 1/2, then \i = 1/2, ( pLOGj ) and ( ]107j ) reduce, 
respectively, to 



da(t) = -iuja(t)dt - na\t)dt + dW t , 
da\t) = ibJo)(t)dt - na(t)dt + dW t . 



(109) 
(110) 



Although dWt and dWt are commutative, we have the 
conservation of the equal-time canonical commutation re- 
lation 



Applying the bra- vacuum ((1| to (|10q ) and (107), we 
have, for any value of u, 

d((l\a(t) = -iw((l\a(t)dt 

-K((l\a(t)dt + V2K((l\dB t , (112) 
d((l\a!(t) = iuj((l\a){t)dt 

-n((l\a\t)dt + V2^((l\dB\. (113) 

Note that these Langevin equations for the vector ((l\a(t) 
and ((l|a^(t) have, respectively, the same structure as (g) 
and §. 



C. Unitary Time- Evolution 

The unitary martingale operator satisfying (|K 
^iven by 



: dMi 



iV2n : (a)dB t - dBja + t.c. 
i (y*dW t + ^dW t 



(114) 



with t.c indicating tilde conjugate. Note that there is no 
cross term between tilde and non-tilde operators. Here, 
we introduced new random force operators 



dW< 



Ik (dB\ 



dBt 



(115) 



which annihilates the bra- vacuum (| of the irrelevant sys- 
tem: 



(\dW? = 0, (\dW? = 0, 
and satisfies the commutation relation 
[dW u dWf] = 2k. 



(116) 



(117) 



The expression (114) is consistent with the microscopic 
Hamiltonian of the linear dissipative coupling. 

The martingale operator (114) satisfies the fluctuation- 
dissipation theorem (p5|) with (|90|). Therefore, we con- 
clude that there exists the system of stochastic differen- 
tial equations with unitary time-evolution operator con- 
sistent with the quantum master equation (|2^) with the 
hat-Hamiltonian (pq) . 

The Langevin equations for a(t) and a'(t) are given by 



da(t) = -iwa(t)dt + V2k dB(t), 
da\t) (t)dt + V2~KdB^(t), 



(118) 
(119) 



where the operators dB(t) and dB^{t) in the Heisenberg 
representation are defined by 
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2k dB(t) = Uj 1 ^) o dB t o U f (t) 

= V2k dB t - m(t)dt, (120) 
2k dB\t) = Uj 1 ^) o dB\ o [>/(*) 

= V2k~ dBj - Ka\t)dt. (121) 



In deriving (120) and (121), we used the properties 



[U f {t) ° V2~K~dB t } = [U f (t), V2k~ dB t ] 



+ -[dUf(t), V2k dB t ] 
KaUf(t)dt, 



(122) 



n R = —i -k (x - x) (p + p) , 

1 - 2 

IJd = — — Kmui(l + 2n) (x — x) 



l\0(t))=-iH\0(t)), 



(131) 



Here, we neglected the diffusion in x-space. The 
Schrodinger equation 



(132) 



gives the quantum Kramers equation p6| . 

The Heisenberg equation for the dissipative system is 
given by 



and 



[U f (t) ° V2k~ dBj] = Ka^Uf(t)dt. 



[Uf{t), dB t ] = [U f (t), dBj] = 0, 



(123) 



(124) 



which comes from the characteristics of the Ito multipli- 
cation: 



(\Uf(t)dB t \) = (\U f (t)dBl\)=0. 



(125) 



Substituting (|12Gj) and ( pL 2 1|) , and applying (|llg| ) 
and (119) reduce to (|112| ) and (113) having the same 
structures as (||) and (Q), respectively. 



-x(t) = i[H(t), x(t)} 

—p(t) = -muj 2 x(t) ~ i/s{p(f) +p(t)} 



(133) 



+iKmuj(l + 2n) {x(t) - x(t)} . (134) 

Applying the bra- vacuum (1| of the relevant system, 
we have the equations for the vectors: 

— <l|p(t) = -mu 2 {l\x{t) - K(l\p(t). (135) 



V. APPLICATION TO QUANTUM KRAMERS 
EQUATION 

A. Quantum Master Equation 

Let us find out the general structure of hat- 
Hamiltonian which is bilinear in (x, p, x, p). x and 
p satisfies the canonical commutation relation 



[x, p] = i. 
Accordingly, x and p satisfies 

[x, p] = -i. 



(126) 



(127) 



The conditions, (iH)~ = iH, and {1\H = give us the 
general expression 



H = H S + in, 



(128) 



where 



~ 1 711 CjJ '~ 

H S = H S -H S , H s = —p 2 + —x 2 , (129) 

Zm 2 



n = n R + n D , 



(130) 



B. Non-Unitary Time- Evolution 

The stochastic Liouvillc equation within the Ito calcu- 
lus becomes 



d\Q f (t))) = -m u dt\o f (t))), 

with the stochastic hat-Hamiltonian 
H ftt dt = Hdt + dM t . 



(136) 



(137) 



Here, the martingale operator dM t satisfying fl3q ) is de- 
fined by 



with 



dM t 



dX t = 



(x - x) (dX t + dX^j 



y/KmUJ 



dB t + dB}) , 



(138) 



(139) 



where dB t , dBj and their tilde conjugates are the oper- 
ators representing quantum Brownian motion (see Ap- 
pendix Q). The generalized fluctuation-dissipation theo- 
rem is given by 



with 



dM t dM t = -2II D dt. 



(140) 
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Taki ng a random average, the stochastic Liouville equa- 
tion ( |l36|) reduces to the quantum master equation (|132| ) 
with |0(f)) = (|0 ; (t)». 

The stochastic Heisenberg equation (the Langevin 
equation) for this hat-Hamiltonian is given by 

dx{t) = i[H f (t)dt, x{t)]-dM{t) [dM(t), x(t)] 

= —p(t)dt + ]-K{x{t) - x(t)}dt, (141) 

TO 2 

dp(t) = -muj 2 x{t)dt - -k {p(t) + p(t)} dt 



(142) 



(143) 



- [dX t + dX t 
where we used the properties 

dX(t)=dX t , dX(t)=dX t . 



Applying the bra vacuum ((1| to (141) and ( |142| ), we 
have the Langevin equations for vectors 



d((l\x(t) = -((l\p(t)dt, 

TO 

d((l\p{t) = -mu? '■((l\x(t)dt - n((l\p(t)dt 
-2((l\dX t . 



(144) 



The aver aged equa tion of motion is given by applying 
|0) to (|144|) and fll45|) in the forms 



|«p(*)» 



(146) 
(147) 



where ((• • •)) = (1|(| • ■ ■ |)|1). The vacuums (| and |) are 
introduced in Appendix O. T hese aver aged equations 
can be also derived from ( p~33| ) and (134) by taking the 
average ((•■■)). 



C. Unitary Time- Evolution 

The martingale operator representing position-position 
interaction may be given by 



dM? = xdX t - xdX t . 



(148) 



We did not include the crossing terms between tilde and 
non-tilde operators to be consistent with the microscopic 
interaction Hamiltonian. 

The fluctuation-dissipation theorem for this martin- 
gale operator is given by 



dMfdM^ = -2II u dt, 



with 



. — (l + 2n) (x-xf. 



(149) 



(150) 



7i u ft dt = H u dt 



dM+ 



where 



h s + in L 



(151) 

(152) 
The mas- 



is the hat-Hamiltonian for the master equation, 
ter equation is different from (p3|). 

The stochastic Heisenberg equations (the Langevin 
equations) for x(t) and p(t) become 



dx(t) = —p(t)dt, 
m 

dp(t) = —mu! 2 x(t)dt — dX t , 
where we used the fact 

dX(t) = dX t . 



Applying ((1| to ( Jl 53| ) and (154), we have the 
equations for the vectors ((l\x(t) and ((l\p(t) in 

d«l|a:(*) = -((l\p(t)dt, 

TO 

d«l|p(*) = ~mu 2 ((l\x{t)dt - ((l\X t , 



(153) 
(154) 

(155) 

Langevin 
the forms 

(156) 
(157) 



(145) which are different from (144) and fl!4E|) 



VI. SUMMARY AND DISCUSSION 

Within the system of non-unitary time-evolution gen- 
erator (non- unitary system), the time-evolution gener- 
ator Vf(t) is constituted by the commutative random 
force operators dWt and dWt- Therefore, the random 
force operators dW(t) and dXif) in the Heisenberg rep- 
resentation is, respectively, equal to dWt and dX t in the 
Schrodinger representation, i.e., 



dW(t) = dW t , dX(t) = dX t . 



(158) 



In the application of the system of unitary time- 
evolution generator (unitary system) to the damped har- 
monic oscillator where the martingale operator is consti- 
tuted by non-commutative random force operators man- 
ifesting the linear dissipative coupling between the rele- 
vant and irrelevant sub-systems, the random force opera- 
tors in the Heisenberg representation are related to those 
in the Schrodinger representation by 



dW{t) = dW t - n lv {t)dt, 
dW l {t) = dWf - K 7 *(t)dt. 



(159) 
(160) 



Then, the Ito stochastic hat-Hamiltonian becomes 



The second terms show up because of the non- 
commutativity. The appearance of these terms is essen- 
tial in order to make the unitary system consistent with 
corresponding master equation. 

On the contrary, in the application of the unitary sys- 
tem to the quantum Kramers equation where the mar- 
tingale operator is constituted only by commutative ran- 
dom force operators manifesting the position-position 
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coupling between the relevant and irrelevant sub-systems, 
the random force operators in the Heisenberg representa- 
tion is equal to those in the Schrodinger representation, 



i.e. 



dX(t) = dX t 



(161) 



Therefore, the unitary system cannot be consistent with 
corresponding master equation. 

The above applications tells us that the origin of dis- 
sipation cannot be quantum mechanical. In spite of this 
unsatisfactory nature of the unitary system, it is attrac- 
tive since hat-Hamiltonian for microscopic system is Her- 
mitian and there is no mixing terms between tilde and 
non-tilde operators. The hat-Hamiltonian should have 
the structure 



H = H — H. 



H ] = H, 



(162) 



for microscopic systems. In fact, we succeeded to extract 
the correct stochastic hat-Hamiltonian for the stochas- 
tic Kramers equation by an appropriate coarse graining 
of operators (the stochastic mapping) in time and corre- 
sponding renormalization of physical quantities f3§f| . The 
simple limit |^9| does not give us the correct Kramers 
equation. This something touchy situation should be in- 
vestigated based on the unified system of stochastic dif- 
ferential equations shown in this paper. It will be re- 
ported in the future publications. 
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APPENDIX A: NEW ASPECT FOR A DAMPED 
HARMONIC OSCILLATOR 

1. Thermal Doublet 

Let us introducing the thermal doublet notation by 

a{ty =1 = a(t), a{ty= 2 = o^t), (Al) 



a{t) 



a T (i), a(t) 



^=2 _ 



-S(t). 



(A2) 



Then, the canonical commutation relation can be written 
as 

[o(t) /i , a(t) v ] = 6"". (A3) 

Note that 

a(ty = V-^t^Vit), d(ty = V- l (t)a^V{t). (A4) 

Making use of the thermal doublet notation, the hat- 
Hamiltonian (p^) reduces to 



H = ujd^a^ +ill + u, 
n = -naPA^a" + k. 



with 



l + 2h -2n 
2(1 + ft) -(l + 2n) J ■ 



(A5) 
(A6) 



(A7) 



The Heisenberg equations for the semi-free particle be- 
come 



-a(t)»=i[H(t), a(ty] 



-i [uS"" - inA^\ a(t) v , 



-a(ty=i[H(t), a(ty] 



a(t) v i [oj6^ - inA^} 



(A8) 



2. Annihilation and Creation Operators 

Let us introduce the annihilation and creation opera- 
tors, 



7 (^= 1 = 7 (i), 7 (*)"= 2 =7*(t), 



(A9) 
(A10) 



by 



1 (ty = B(t)^a(t) u , 7(t)" = Sp- 1 !*)"', (All) 
with the time- dependent Bogoliubov transformation: 



B(tr = [ 1+ ^ (<) 



(A12) 



where n(t) is the one-particle distribution function satis- 
fying the Boltzmann equation (87). 

The annihilation and creation operators satisfy the 
canonical commutation relation 



[y(tr, 7(*)1=<^, 



(A13) 



and annihilate the bra- and ket-vacuums at the initial 
time: 



7 (t)|0) =0, (l|7*(t) = 0. 



(A14) 



The equation of motion for the thermal doublet 7(t) M 
is derived as 



^l a (tr + B(tr^ a (tr 

dB(t) , 



dt 



B-'it) 



■i [B{t) {lu l~~iKA)B- 1 {t)Y u 1 {tY 



= -i [uS^ - iKT^ v ] 7(t)", (A15) 
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where the matrix 73 is defined by 

^22 -, 12 21 n 

For the third equality, we used the relations 
dB{t) ( 1 -1 A dn(t) 



dt 
dB{t) 







dt 



S(t)AB- x (t) = r 3 + 2 [n(t) - ft] r+, 



where 



1 




(A16) 

(A17) 

(A18) 
(A19) 

(A20) 



The Boltzmann equation ( CT ) has been used also. 
The solution of (A15) is given by 

7(t)" = exp { — z (w<5^ - iivrf) (t - t')} >y(l?) v . (A21) 

This expression gives 



(l\~f(thHt')\0) 



— i{uj — lK){t — t') ^ 

— iK)(t— t ) 

1|7(*')7*(*)|0) = <l| 7 (f)7*(*)|0) 

_ —i(ui+iK)(t—t') 



l| 7 (f)7*(f)|0) 



(A22) 
(A23) 



3. Two-Point Function (Propagator) 

The time-ordered two-point function G(t, t'Y u has the 
form 



G{t 1 t'Y u = -i{l\T[a{tyd(t'Y] |0) 

= [B-\t)g{t,t')B{t')Y v , 



(A24) 



where 



G(t,t'Y" 



-i(l\T tyWW] 10) 
(G R (t,t>) 
1. G A (t,t>) 



with 



G R (t,t') = -i9(t - t^e-^-^-t"* , 
G A (t,t') = i6{t' - t)e- l ^+ lK ^ t '\ 



(A25) 

(A26) 
(A27) 

In deriving the above expression, we used the elements 
of the solution (A.21) with some algebraic manipulations. 
For example, 

g(t,t') n = -i(i\Th(th(t'f]\o) 

= -i{(5>(t-O<l|7(*)7*(*')|0) 

+e(t'~t){i\ 1 Ht'h(t)\o)} 

= -i6{t-t')e- t{uJ - lfi){t - t ' ) 
= G R (t,t'). 

In the third equality, we used ( |A2^) and ( jlgf ). 



(A28) 



4. Miscellaneous 

The representation space (the thermal space) of 
NETFD is the vector space spanned by the set of bra 
and ket state vectors which are generated, respectively, 
by cyclic operations of the annihilation operators 7(f) 
and 7(4) on (1|, and of the creation operators 7+(t) and 
7*(t) on |0>. 

The normal product is defined by means of the annihi- 
lation and the creation operators, i.e. 7*(t), 7 stand 
to the left of 7(f), 7(f). The process, rewriting physical 
operators in terms of the annihilation and creation oper- 
ators, leads to a Wick-type formula, which in turn leads 
to Feynman-type diagrams for multi-point functions in 
the renormalized interaction representation. The inter- 
nal line in the Feyn man-t ype diagrams is the unperturbed 
two-point function ( A24 ). 



5. Condensation of Particle Pairs 

Introducing the annihilation and creation operators in 
the Schrodinger representation 



r =1 =7t, 



7^ = 7* 



(A29) 
(A30) 



7 = -It, 
by the relation 

7 (tr = v~\Mv{t), j(tY = v-\t)^v(t), 

(A31) 



with V(t) being specified by (p9[), we can rewrite the hat- 
Hamiltonian (pq) as 



H = lj (7*7* - 7*7t ) - ill, 



(A32) 



with 

n 



-K (7*7* + 7% + 2 [n{t) - n] 7*7*) . (A33) 



It is easily deriv ed by means of the doublet notation 1 

Substituting ( A32| ) into the quantum master equation 
(pq), we have 



d_ 

at 



|0(t)) =-2«[n(t)-n]7*7*|0(t)) 



dn{t) _ 
dt 



7*7*|0(*)}. 



(A34) 



It is solved to give 
|0(t)) = exp 



* ,dn(f) 



—7 7 



|0) 



= exp [n(t) - n(0)] 7*7* |0) 



(A35) 
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This expression tells us that the time-evolution of the un- 

o o 

stable vacuum is realized by the condensation of j^jZ- 
pairs into the vacuum. The attractive expression ( |A35| ), 
which was obtained first in p9fl , led us to the notion of 
a mechanism named the spontaneous creation of dissipa- 
tion 0,|5U|^^] . The corresponding order parameter is 
given by 



(lfofc7t|0>=-n(t)+n(0) 



where we used the relation 



7t = 7t=o ~ [n(t) - n(0)]7 



(A36) 



(A37) 



We can obtain the results (A35) and ( A36| ) only by al- 
gebraic manipulations. This technical convenience of the 
operator algebra in NETFD, which is very much similar 
to that of the usual quantum field theory, enables us to 
treat open systems in far-from-equilibrium state simpler 
and more transparent [f30| -[35[. 

It also shows that the vacuum is the functional of the 
one-particle distribution function n k (t). The dependence 
of the thermal vacuum on Ufc(t) is given by 



|0(*))=7^|0(t)) 



Sn k (t) 



(A38) 



We see that the vacuum |0(t)) represents the state where 
exists the macroscopic object described by the one- 
particle distribution function n k (t) . The master equation 
(p3|) can be rewritten as 



d_ 

at 



d^k 



dn k (t) 



dt Sn k (t) 



\0(t)) = 0. (A39) 



This shows that the reference vacuum, in this case, is 
migrating in the super-representation space spanned by 
the one-particle distribution function {n k {t)} with the 
velocity {dn k (t) / dt} as a conserved quantity. 



It is easy to see from the normal product form (A32) 
of H that it satisfies (|l^) , since the annihilation and cre- 
ation operators satisfy 



7t|0(*)> =0, (1|7* = 0. 



(A40) 



The hat-Hamiltonian 
form 



y can be also written in the 



H = u (dU - cftd) - in (dU + J f d) , 



(A41) 



where d» =1 = d, eP= 2 = <$ and d» =1 = d+, d»= 2 = -d 
are defined by 



d" = B^a v , # = a v B- lv »< 



with 



B^ 



1 + n —n 
-1 1 



(A42) 



(A43) 



The ket-thermal vacuum, |0) = |0(0)), is specified by 
@ which can be expr essed in terms of d and cT, which 
are introduced in ( |A42 ) below, as 



d|0) = (n-n) d f |0). 



(A44) 



It is easy to see from the diagonalized form (A41) of 
H that 

d(t) = V^it) d V(t) = d e -(^+ K )*, (A45) 
J tt (i) = V^it) $ V(t) = i? e-^-^K (A46) 

The difference between the operators which diagonalizc 
H and the ones which make H in the form of normal 
product is one of the features of NETFD, and shows the 
point that the formalism is quite different from usual 
quantum mechanics and quantum field theory. This is 
a manifestation of the fact that the hat-Hamiltonian is 
a time-evolution generator for irreversible processes. In 
thermal equilibrium state, n(t) — n, they coincide. 



6. Irreversibility 

Let us check here the irreversibility of the system. The 
entropy of the system is given by 

S(t) = -{n(t)ln.n(t) - [l + n(t)]ln[l + n(t)]}, (A47) 

whereas the heat change of the system is given by 

d'Q = udn. (A48) 

Thermodynamics tells us that 

dS = dS P . + dSi. dS,, = d'QjT n - (A49) 



dSi, dS e = d'Q/T R , 
dSi > 0. 



(A50) 



The latter inequality ( A50| ) is t he se cond l aw of thermo- 
dynamics. Putting ( |A47| ) and flA48| ) into QA49p , for dS 
and dS e , respectively, we have a relation for the entropy 
production rate [|53l] 



dt 



dS 
~dt 



dS e n(t)[l + n] 
—— = 2k \n(t) — n\ in — —-^ 

dt 1 w J n[l + n(t)} 



> 0. 



(A51) 



It is easy to check that the expression on the right-hand 
side of the second equal ity s atisfies the last inequality 
which is consistent with ( A5C ) . The equality realizes ei- 
ther for the thermal equilibrium state, n(t) — n, or for 
the quasi-stationary process, k — > 0. 



APPENDIX B: ITO AND STRATONOVICH 
MULTIPLICATIONS 

The definitions of the Ito jl7| and the Stratonovich | 
multiplications are given, respectively, by 
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X^(t)-dY^(t) 
= X iH \t) \Y<- H \t + dt) - Y {H) (t) 

dxW(t)-yW(t) 

= xW(t + dt)-XW(tj\YW(t), 



(Bl) 
(B2) 



and 



Y^{t + dt) -Y {H \t) 
dX^ (t) o Y {H) (t) = \x^ H 3 (i + dt) - (i) 

r( H )(t + dt) + y( H )(t) 



(B3) 



(B4) 



for arbitrary stochastic operators X^ H \t) and Y^^t) 
in the Heisenberg representation. From (|Bl|), (B2) and 
(B3), (B4), we have the formulae which connect the Ito 
and the Stratonovich products in the differential form 

i< ff )(t) o dY^(t) = x( ff )(i)dy( ff )(i) 

+-dx^(t)-dy^(i), (B5) 
diW(t) o yW(t) = dx(*)(t) ■ rW(t) 

(B6) 



The connection formulae for the stochastic operators 
in the Schrodinger representation are given, in the same 
form as (B5) and (pq), by 



+idx( s )(t)-dy( s )(0, 
+idx( s )(t)-dy( s )(t), 



(B7) 
(B8) 



where the operators X( s \t) and dX( s '(t) in the 
Schrodinger representation are introduced respectively 
through X( H \t) = V^tyX^^Vfit) and dX^t) = 

v f -\t)dx^(t)v f (t). 



APPENDIX C: QUANTUM BROWNIAN 
MOTION 

Let us introduce the annihilation and creation opera- 
tors bt, b\ and their tilde conjugates satisfying the canon- 
ical commutation relation: 

[h, b\,] = 5(t - t'), [bt, b\,] = S(t - U). (CI) 

The vacuums (| and |) are defined by 



6t|) = 0, &t|) = 0, (\b] = (\b t . 

The argument t represents time. 
Introducing the operators 

rt—dt 



(C2) 



pz—az pi 

B t = dB t > = / dt' b v , 

Jo Jo 



BJ 



t-dt 



dBj 



dt' b\„ 



(C3) 
(C4) 



and their tilde conjugates for t > 0, we see that they 
satisfy B(0) = 0, B*(0) = 0, 



[B s , B\] = min(s,t), 



(C5) 



and their tilde conjugates, and that they annihilate the 
vacuum |) with the thermal state condition for (|: 

dB t |) = 0, dB t |) = 0, {\dB\ = {\dB t . (C6) 

These operators represent the quantum Brownian mo- 
tion. 

Let us introduce a set of new operators by the relation 
dC? = B^ v dB v t , (C7) 
with the Bogoliubov transformation defined by 



1 + n —n 
-1 1 



(C8) 



where n is the Planck distribution function. We intro- 
duced the thermal doublet: 

dB? =1 = dB t , dB?= 2 = dfy, (C9) 
dB? =1 = dB\, dB?= 2 = -dB t , (CIO) 

and the similar doublet notations for dC^ and dC? . The 
new operators annihilate the new vacuum (|, and have 
the thermal state condition for |): 

dC t \} = 0, dC t \) = 0, {\dC\ = {\dCt. (Cll) 

We will use the representation space constructed on 
the vacuums (| and |). Then, we have, for example, 

(\dB t \) = (\dBt\) = 0, (C12) 
(\dB\dB t \) = ndt, (\dB t dB\\) = (1 + ft) dt. (C13) 



They can be written 

dB^ t dB t = dB t dB t = ndt, 
dB t dBj = dB\dB\ = (1 + n)dt, 

as weak relations. 



(C14) 
(C15) 
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